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A rigorous theoretical framework within which quarkoniuni production can 
be studied is provided by the NRQCD factorization formalisniEl. This approach 
is based on NRQCD, an effective field theory that that can be made equivalent 
to full QCD to any desired order in the heavy-quark velocity v. The NRQCD 
Lagrangian is 
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where is the gluon field-strength tensor, is the Dirac spinor field for 
a light quark, = + igA^ is the gauge-covariant derivative (with g be- 
ing the QCD coupling constant), toq is the heavy quark mass, and tp and x 
are 2-component fields that annihilate heavy quarks and create heavy anti- 
quarks respectively. Color and spin indices on the fields have been suppressed. 
Eq. |l| describes a fully relativistic field theory for the light degrees of freedom 
coupled with a Schrodinger field theory for the heavy quarks and antiquarks. 
The relativistic effects of full QCD are reproduced through additional terms 
represented by the • • • in Eq. |l|. In principle there are infinitely many terms. 
However, using NRQCD w-scaling rules it is possible to retain only those that 
contribute through a given order in v. 

The NRQCD factorization formalism is a rigorous derivation within NRQCD 
of a factored form for quarkonium production and decay rates. A central result 
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is that inclusive quarkonium production cross sections must have the form 



a(A + B - H + X) = £ (°n) , (2) 

n m Q 

where the represent NRQCD four-fermion production operators, with 
the index n labeling color and angular-momentum quantum numbers. The 
NRQCD matrix elements (O^) = (0|O^|0), parameterize the hadronization 
of a heavy-quark-antiquark pair with quantum numbers n into a quarkonium 
state H . They scale as the parameter v raised to a power. The short-distance 
coefficients F n are obtainable in perturbation theory as a series in a s (mQ). 
Information on the order in v of the matrix elements, along with the depen- 
dence of the F n on coupling constants, permits us to decide which terms must 
be retained in expressions for observables to reach a given level of accuracy. 

Consider the specific case of J r /ip production. In many instances the most 
important NRQCD matrix elements are (Of ( 3 Si)), (Of ( 3 Si)), (O'gfSo)}, and 
(Of ( 3 Pj)) ■ Inclusive production of J/ip from &-decay provides two measurable 
combinations of these matrix elements. The first one is the helicity-summed 
rate T(b — » J/ip + X). The second combination concerns the polarization 
parameter a appearing in the electromagnetic decay rate of J/ip to lepton 
pairs: 

-U^fi+n-(8)) cxl + acos 2 6l, (3) 

a cos 

where the polar angle 9 is defined in the J/ip rest frame for which the z-axis 
is aligned with the direction of motion of the J/ip in the lab. 

The branching ratio Jor b — > J/ip + X is calculated using a matching 
procedure El. One obtains U 

BR(b — » J/ip + X) = 

O.OO2(0f ( 3 ^)) +O.1(0^( 3 5 1 )) +O.2(0^( 1 5o)) +0.6 ^^ (3 f°^ .(4) 

Only the leading color-singlet piece and the leading color-octet pieces in the 
relativistic i>1ij-expansion are considered. The above formula concurs with pre- 
vious results □. 

According to the NRQCD w-scaling rules, the color-octet matrix elements 
in Eq. || are all expected to be suppressed by v 4 with respect to (Of ( 3 <Si))- 
However the short-distance coefficients in the color-octet terms are some 50 
times larger than the color-singlet coefficient! 
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Beneke and RothstcinH and Braaten and ChenB have developed techniques 
for deriving the production rates of quarkonia with specified helicities. Apply- 
ing their methods one obtains an expression for the polarization parameter! 3 ] 

-0.39(Of ( 3 ^)) - 17(Ot( 3 S 1 ))+52(Ot( 3 P ))/ml 
(Of (350) + 44(0^(3^)) + 6l(O?(iS )> + 211(0^ (3p ))/m2 ' 

While experimental determinations of helicity-summed BR(b — * J/tp + X) 
have already been carried outQ, a measurement of the polarization parameter 
a is not yet available. Anticipating the availability of this latter measurement, 
it is interesting to determine the range of a which is consistent with existing 
information on the matrix elements (Of ( 3 Si)), (Of( 3 Si)), (O^( x 5 )>, and 
(Og ( 3 Po)) r Zm^, and with various constraints on linear combinations of these 
quantities u. The maximum value for a is —0.09, and the minimum value for 
a is —0.33. 

Unfortunately, due to the large number of poorly determined parameters 
the expected range of a as predicted by NRQCD is large. It may be that 
an accurate NRQCD prediction of the polarization parameter (beyond leading 
order) will not be possible for a long time. This possibility, however, does 
not degrade the importance of this calculation, since an experimental deter- 
mination of a — in conjunction with these results — will most certainly serve 
to tighten the constraints on the possible values of the color-octet matrix el- 
ements. In fact, the NRQCD prediction for a is very sensitive to the values 
of the color-octet matrix elements. This offers the hope that a measurement 
of the polarization of J/ tp produced in b — > J /if) + X will be instrumental in 
determining ((^( 3 Si)>, (O^So)), and (Ot( 3 P ))- 
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